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Abstract 

I propose that the properties of QCD perturbation theory should 
be investigated when the boundary state ('perturbative vacuum') at 
t = ±oo includes gluons. Any boundary state that has an overlap with 
the true QCD ground state generates a perturbative series that (when 
summed to all orders) is formally exact. Through an analogy with the 
boundary condition corresponding to a fermion condensate, I propose 
an explicit form for a 'perturbative gluon condensate' that suppresses 
low momentum gluon production, thus generating an effective mass 
gap. Standard perturbative calculations are modified only through a 
change in the ie prescription of low momentum (\p\ ;$ Aqcd) gluon 
propagators. Gauge invariance is expected to be preserved since this 
modification is equivalent to adding on-shell external particles. Renor- 
malizability is unaffected since only low-momentum propagators are 
modified. Due to the asymptotic low momentum gluons boost invari- 
ance is not explicit. Lorentz invariance should be restored in the sum 
to all orders in analogy to standard bound state calculations. 
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1. Introduction 

Hadron wave functions appear phenomenologically to be frame dependent. 
In the infinite-momentum (or light-cone) frame the proton is observed to 
have, in addition to its uud valence quarks, important gluon and sea quark 
components jl]. In particular, gluons carry about half of the proton momen- 
tum. This measurement of the proton is rather rigorously justified by the 
QCD factorization theorem ||. The success of perturbative QCD predictions 
for a large number of hard scattering processes has established QCD as the 
correct theory for the strong interactions. 

The non-relativistic quark model (NRQM) p], |3], f|] provides a less rigor- 
ous but phenomenologically very successful and simple rest frame picture of 
hadrons as non-relativistic bound states of 'constituent' quarks. The masses 
of the u,d constituent quarks are (9(300 MeV), which is considerably larger 
than the 'current' quark masses which are relevant for short-distance pro- 
cesses. The 'missing' gluon and sea quark degrees of freedom appear to be 
frozen in the structure of the constituent quarks. 

The simple regularities of the hadron spectrum, coupled with the success 
of QCD as applied to hard processes, obviously invites efforts to find a QCD 
justification of the NRQM (see, eg, [[| |, [], [7|, The challenge is to find a 
formulation which in a first approximation retains the simplicity of the quark 
model, yet allows corrections specified by QCD to be evaluated to arbitrary 
order. 

The similarities of the hadron spectrum with QED bound states, together 
with the success of perturbative calculations in that theory, suggests the use 
of a perturbation expansion also in QCD. Such an expansion is determined by 
the lagrangian and by the boundary conditions at asymptotic times. Central 
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properties of QCD like the finite range of the color force and confinement are 
commonly associated with a non-trivial ground state of the theory, the 'gluon 
condensate' |@, [K| . In this paper I shall propose a specific boundary condition 
on QCD perturbation theory which is motivated by the gluon condensate and 
by a suppression of soft gluon production. The purpose is not to model the 
condensate in detail. Just as in QED it may suffice to achieve so much overlap 
with the true ground state that low orders of the perturbative series already 
incorporates the main physical features of the theory. Corrections are then 
given systematically by the higher orders of the expansion. 

This work points to a perturbative expansion of QCD which appears not 
substantially more difficult to evaluate than the standard one, but which has 
a number of novel features. The usefulness of this approach can only be 
judged after a further study of the properties of that expansion. 

2. A fermion condensate 

The asymptotic states that we impose on perturbative expansions at initial 
and final times (t = ±00) should have an overlap with the true ground 
state of the theory. This guarantees that the full perturbative expansion 
formally gives exact results. In euclidean formulations this fact is particularly 
clear since the time development of energy eigenstates is given by exp(— Er), 
implying a dominance of the true ground state (of lowest energy E) in the 
limit t — > 00. In minkowski space the same result is obtained using an is 
prescription^ 

This freedom in the choice of boundary states allows for a whole set of for- 
mally equivalent perturbative expansions. Since all expansions are expected 
to diverge, their equivalence is of more formal than practical significance. 
2 For a discussion of boundary states in field theory see, eg, Ref. [O 
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From a practical point of view what matters is that the lowest orders already 
incorporate the main physical characteristics of the theory. 

We have little understanding of the structure of the QCD gluon conden- 
sate in terms of Fock state wave functions. It nevertheless seems plausible 
that the vacuum wave function components involving gluons and quarks of 
3-momenta smaller than the characteristic QCD scale Aqcd are strongly 
modified. In this respect, the gluon condensate may resemble a fermi con- 
densate with fermi momentum of O(Aqcd)- In a fermi condensate the ex- 
clusion principle prevents pair production below the fermi momentum. It 
seems desirable to have a similar property for gluons, to suppress soft gluon 
production which can give rise to long-range color correlations. 

In this section I recall how perturbation theory is modified in the presence 
of a fermion condensate. Only the is prescription is affected - which is enough 
to have significant consequences. I shall then in the next section use this 
as a guide for constructing a 'perturbative gluon condensate', namely one 
that results in an analogous modification of the is prescription for the gluon 
propagator. Having shown that there exists a boundary condition which 
implies such an is modification for gluons it can for many practical purposes 
be forgotten, and the usual feynman diagrams be evaluated with modified 
(low momentum) propagators. 

The standard free fermion propagator 




(1) 



is in momentum space 



Sf(p) 





(2) 



> 2 — m 2 + ie 



(p°-E p + ie)(p° + E p -ie) 
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If we add an antifermion to the initial and final states, 



{Q\dy(k')T^(x)^(y)}d\(k)\0) = tS F (x-y)2E k (2n) 3 5 3 (k-k>)5 xx , 

+ v(X', k')v(X, k)e ik ' ' x ~ ik ' v , (3) 

the feynman propagator is multiplied by the annihilation amplitude for the 
inserted antifermions, and there is a new term corresponding to a mixing of 
the antifermion propagating from x to y with the antifermion in the in- and 
out-states. 

For a condensate we would fill both helicity states at a given momentum 
k. The free propagator 

iS(x -y) = (0|d 1/a (fe)d_ 1/2 (fe)T[V(x)^(y)]dL 1/2 (^)4/2(^)|0) ( 4 ) 
is then in momentum space 

S(p) = [(2.) W(0)] 2 { :| (5) 

where 

a ( \ = t +m (R\ 

E{P> ( p °-E k + ie)( P ° + E k + ie) U 
differs from the feynman propagator only in the ie prescription at p° = —E k . 
Since the antifermions inserted in the definition (f|) are on-shell, it is clear 
that a mixing between them and the propagating fermion only can occur at 
the antifermion pole of S(p). Adding antifermions for all momenta \k\ < A, 
the corresponding propagator will equal Se(p) for all \p\ < A. 

The addition of (anti)fermions at t — ±oo influences the fermion propa- 
gators in feynman diagrams at any order of perturbation theory exactly as 
it does the lowest order propagator above. This can be easily seen using the 
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generating functional of green functions (in a theory like QED or QCD), 



z[C,C; J] = ex P 



iS{ 



hit 



'55 5 ' 



z B [j]z F \U] 



(7) 



where Si n t is the interaction part of the action and Zb, Zp are free functional 
of the boson (J) and fermion ((, () sources, respectively. The free fermion 
functional is 

d 4 p - 



z f[CX] = exp 



(2vr) 



-X(-p)s F (p)C( P ) 



(8) 



with the feynman propagator given by Eq. (fj). 

It is instructive first to rederive the result ([5]) for the free propagator with 
the boundary states @ using the generating functional. The propagator 
Sf(p) is diagonalized by the sources z, z of definite helicity (A) and energy 
signature (±), 

, o 
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2E„ A 



[u{\,p)z\(p)+v(\-p) z x _(p) 



C(p) = 

V 

a-p) = iE^H^M + ^-M^H 

In the new basis we have 



(9) 



ZAz,z)=exv\ij^Y: 



(2tt)< 



z +{~P) z +(p) _j_ z -{—p)Z-(p) 



(10) 



p° - E p + is p° + E p - is 

The generating functional Ze of the modified propagator (§) (for some given 
k) differs by the sign of ie in the second term of Eq. (|I0"D. Hence (I suppress 
the 3- momentum p and factors (27r) 3 2i?5 3 (0) in the following), 

ex P i I ■7t-C(-p)Se(p)((p) 

U[l + z x _(E)z x _(-E)]z F [z,z] , (11) 



Z E [z,z] 
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where I used (p° + E + ie)^ 1 — (p° + E — 
exp(^2:) = 1 + zz for grassmann sources z, z. 
As a function of time, 



IE 



,-1 



-2m5(p° + E), and 



z(p°) = 
z(-p°) = 

the free generating functionals are 



dt' z{t')e it,p ° 
dt" 'z(f)e- it " p ° 



(12) 



Z F = expl^dt'dt"J2[zi(t"W -t')e 

-z x _(t")6(t' - Oe~ iB(t '~ t " ) ^(0] } 



(13) 
(14) 



This expression for Ze can now be compared with the one obtained by ex- 
plicitly differentiating Ze wrt. its sources at t = ±oo, corresponding to the 
boundary states of Eq. (|). One readily finds 

5 2 



lim j~J 



ti—> — oo 



JE(tj-ti)_ 



Sz±(ti)6zl(t f )\ 



(15) 



This result extends immediately to the full interactive functional (0), since 
the derivatives in Eq. (|15|) commute through the derivatives in exp(zSi„ t ). 
This means that using the modified fermion propagator Se of Eq. (0) every- 
where in a perturbative calculation of an arbitrary green function (for some 
given 3-momentum k of the propagators) is exactly equivalent to calculat- 
ing the same green function using ordinary feynman propagators but with 
additional incoming and outgoing antifermions as in Eq. (Q). 

The change of is prescription suppresses fermion pair production at the 
corresponding value(s) of k, as required by the exclusion principle. This 



6 



can be seen directly for, eg, a fermion loop correction to a gauge boson 
propagator. The loop gives no contribution at those values of the fermion 
momenta k at which external fermions have been introduced, since the poles 
in the loop momentum p° then are all below the real axis and the p° integral 
may be closed in the upper half plane. 

Since I have shown that the propagator modification is equivalent to 
adding particles at t — ±00 gauge invariance is likely to be preserved. For- 
mally, the ward identities involve inverse propagators, for which the sign of 
ie is irrelevant. 

3. A boson 'condensate' 

Boundary conditions like that of Eq. (^) with (anti)fermions added to the 
in- and out-states are relevant in situations involving fermion condensates, 
but not for typical applications of QCD. The QCD vacuum has zero baryon 
number, and thus no overlap with states having extra (anti) quarks. The 
propagator modification nevertheless seems phenomenologically interesting 
for gluons, since it suggests a 'freezing' of the low momentum gluon d.o.f.'s. 
Effective gluon and constituent quark masses can be generated through loop 
corrections due to the propagator modification, presumably without loss of 
gauge invariance (but with loss of lorentz invariance order by order, see sec- 
tion 4). 

I shall show that there is a boundary condition which implies an analo- 
gous modification of the ie prescription for boson propagators as the one for 
fermions discusses above. Not surprisingly, this 'perturbative boson conden- 
sate' involves an indefinite number of external bosons. For simplicity, I shall 
consider scalar bosons only. The generalization to real (transverse) gluons 
should be straightforward. 



7 



The free boson functional appearing in Eq. (|7]) is (for scalars) 



Z B [J}= exp 



±p 



(16) 



Since we shall be dealing with the free functional (the generalization to the 
interacting one will again be straightforward), it is sufficient to consider a 
single 3-momentum p, and keep only the bose symmetrization over ±p as 
indicated in Eq. (fTS). The feynman propagator is 

Df ^ = p 2 -m 2 + ie = 2E {p°-E + ie ~ p° + E-ie) ^ 

where E = ^Jp 1 + m 2 . A modification of the is prescription at the p° = —E 
pole gives 



De(p) 



W^EnwTEnrr DM+ iE 5<J, ° +E) (18) 



Note that the same generating functional is obtained if the ie prescription is 
changed instead at the p° = +E pole. Thus 

1 



so that 



De(p) = 



dp 



(p° - E - ie)(p° + E - ie) 



De(-p) 



E / ^J(-P)D E (P)J(P) = E / -^J(-p)d e (p)j(p) 



±p 

In (t,p)-space 



dp 



2tt 



J(p°,p)= JdtJ(t,p)e^° 



(19) 



(20) 



(21) 



we have 



Z B [J] = e W ij2^Jdt'dt"J(t",-^) 

x 0{t" - t')e- iE( t"-^ + 6{t' - t")e iE ^"-^] J(t',p)} 



(22) 



The generating functional for the modified scalar propagator flIB|) is then 

dp 



Z E [J] = exp 



exp 



exp 



E 

±p 



-J(-p°,-p)D E (p)J(p°,#) 



(23) 



-Y,^J(E,-P)J(-E,P) 



Z B [J] 



E 4^/ -p)e iE V'-*)j(t,p) 



Z B [J] (24) 



Eq. ( f2~4"|) may be compared with Eq. (|14"D in the fermion case. Due to 
the grassmann algebra, the exponential factor multiplying Zp contains only 
a single power of the fermion sources z, z. In the boson case the factor 
multiplying Zb in Eq. fl2~3| ) contains arbitrary powers of the sources J. It can 
be reproduced only by differentiating Zb[J] an arbitrary number of times, 
corresponding to an indefinite number of incoming and outgoing bosons. 

A single boson of momentum p in the in-state is obtained as 

szjM = m 

ti-- oo 6J(t h -p) V2E 
where we used lim^-^-oo 9(ti — t) =0. Similarly an outgoing boson corre- 
sponds to 



J dte- tE(t ~ u) J(t,p) Z B [J] 



(25) 



SZ B [J] 

tf^+oo5J{tf,p) 



lim 



1 

2E 



dte -iE{t s -t)j^ _^ 



z B [J] 



(26) 



Having a boson both incoming and outgoing is then given by 
lim e iE ^2E - - ^ff^ - = 

ti->— oo oJ(tf,p)oJ(ti,—p) 



tf— >+oc 



1 + — / dt'dt" J(t", -p)e 



iE{t"-t') 



J(t',p) 



Z B [J] 



(27) 



Further differentiation wrt. </(£», — p) and J(tf,p) now operates also on the 
first factor in Eq. (p7|). However, this gives back the factors in Eqs. ( p5|) 



and , respectively. Hence applying the double derivative of Eq. ( p7|) any 
number of times on Zb generates a polynomial factor in xy, where 



V2E 

According to Eq. 



x = -t=fdt"J(t",-p)e* Et '' 

= Jdt'e- iEt 'j(t',f) . (28) 



Z E [J]=exp\-~^xy)z B [J] . (29) 

We need to consider only how to generate the +p term in Eq. (|2"9"| ) through 
repeated differentiation of Zb as in Eq. (|27|). The — p term will then be 
obtained similarly through repeated 5 2 /5J(tf, —p)5J(ti,p) differentiation. 
The function 

( d 2 \ 

f(xy) = exp I ^q^J exp(xy) (30) 

provides an adequate model for the present problem. As seen from Eq. (|22]) , 
Z B is not of the form exp(xy) due to the ^-functions, but as in Eqs. (^5| - 



27]) Zb acts precisely like exp(xy) when differentiated in the limits t; t — > 
— oo, t/ — > +oo. Hence the polynomial in xy generated by the derivatives in 
Eq. ( |3T)D will be the same as that generated from Zb- 

It is straightforward to evaluate f(xy) in Eq. (KOI) by using the identity 



1 \ I r°° /I 



exp ( -x J = . — J duexp -u + ux ) (31) 



to express 

exp(xy) = exp 



0X1) (_I^ exp /_i y2 ) (32) 



10 



as a three-fold integral. The integral resulting from applying the derivatives 
in Eq. fl3"U| ) is gaussian and gives 

exp ( A ^k) eMxy) = ih^i^h) 

= T^A exp 

Requiring that the first exponent in Eq. ([33|) be —xy/2 according to Eq. ( p9|) 
gives the result A = — 1. 

We have thus shown that the generating functional Z%[J] ( p3|) of the 
scalar propagator De{p) (0), which differs from the feynman propagator 
Dp{p) QT7| ) by the sign of is at the p° = —E pole, is equivalent to the standard 
generating functional Zb[J] (0) of feynman propagators differentiated wrt. 
sources at t = ±oo, 

Z E [J] = 4 exp 

As noted in Eq. (P0|), the same result obtains if the is prescription is modified 
at the p° = +E pole instead. 

The source derivatives in (0) commute through the interaction term in 
the definition (0) of the full generating functional of green functions in the 
interacting theory. Hence the above result establishes that a perturbative 
calculation (to arbitrary order) which uses the propagator De{p) (0), with 
its non-standard ie prescription, is equivalent to a standard perturbative 
calculation using feynman propagators in the presence of a 'perturbative 
condensate' of incoming and outgoing particles as specified by the source 
derivatives in Eq. fl34|) . 
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1-A 



xy exp(xy) (33) 



lim e iE{t f- u) 2E 



-oo 
-oo 



SJ(t f ,p)5J(ti, -p) 



Z B [J] (34) 



4. Discussion 

I have argued that it may be useful to consider perturbative expansions 
of QCD using non-trivial boundary conditions at t — ±00, given that the 
ground state of the theory is a gluon condensate. All expansions in which the 
boundary states overlap the true ground state are formally equivalent and a 
priori equally good. 

I investigated a particular case which is the bosonic equivalent of a fermion 
condensate, and in which the is prescription of low momentum boson prop- 
agators is modified. Such a propagator modification corresponds to a super- 
position of standard perturbative calculations where 0, 1, 2, etc. bosons are 
added both to the initial and final state, as expressed by Eq. (0). I have not 
shown that these boundary states have an overlap with the true QCD vac- 
uum (but then, neither do we know that the standard perturbative vacuum 
has such an overlap). 

The relevance of this expansion depends on its theoretical and phenomeno- 
logical viability, which remains to be demonstrated. Gauge invariance is 
among the important properties that should be explicitly verified. 

Since only low-momentum < Aq C d) propagators are modified, the 
successful results of 'hard' QCD processes remain unaltered. In particular, 
the renormalization procedure will not be affected in any way by the modi- 
fications suggested here. 

The most striking difference compared to standard perturbation theory is 
the lack of boost invariance order by order. Contrary to what might first ap- 
pear, this need not signal a breakdown of lorentz symmetry for the full series. 
The true asymptotic degrees of freedom are the hadron bound states, which 
do not occur at any finite order of perturbation theory. Physical symmetry 
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requirements should be imposed only on resummations of the series. 

The subtleties of lorentz invariance in bound state calculations is known 
from QED. As an example |12 |, consider the lippman-schwinger equation 



G T (E) = K(E) + K{E)S{E)G T {E) (35) 

for the (truncated) green function Gt of a 2 — > 2 process with cm. energy 
E. Iterating this equation generates an expansion of Gt in powers of the 
propagator 5* and the kernel K. While the standard perturbative expansion 
in a is unique (up to renormalization conventions) for the green function 
Gt, this is not so for S and K separately. Rather, we can choose the form 
of the propagator S freely, be it of relativistic (dirac) or non-relativistic 
(schrodinger) form. Eq. ( |3~5l) then determines the corresponding perturbative 
expansion of the kernel K. At a pole of the (full) green function of the form 

G{E) = ^ n ^ n — h regular terms (36) 
E — E n 

the lippman-schwinger equation implies a bound state equation of the form 

S-\E n ) = K(E n )^j n . (37) 

For a non-relativistic propagator S this will have the form of a schrodinger 
equation, but it will give exact results provided the full perturbative series 
for the interaction kernel K is used. 

It should furthermore be realized that the transformation properties of 
equal-time bound state wave functions under lorentz boosts is quite non- 
trivial. The requirement that the constituents should be evaluated at equal 
time in all frames is inconsistent with explicit space-time covariance, even 
for non-relativistic QED bound states. 
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A simple example serves to illustrate the novel aspects of the frame de- 
pendence of equal-time wave functions. There is a bound state equation in 
QED 2 for which it is possible to relate explicitly the solutions in different 
lorentz frames, and thus verify that they have the correct transformation 



properties [fL3]j . The wave function of a two fermion bound state is written 

ip(t,xi,x 2 ) = exp(-iEt) exp (ik — %2 ^j x{xi — x 2 ) , (38) 

where X\,x 2 are the positions of the constituents and t their common time. 
Both the bound state energy E and the 2x2 dirac wave function x depend 
on the bound state cm. momentum parameter k. The bound state equation 
for x is 

- id x [a, x(x)} + hi {a, x(x)} + mi7°x(x) - m 2 x{x)rf = {E - V(x)) x(x) 

(39) 

where mi,m 2 are the constituent masses and V(x) = |e 2 |x| is the instan- 
taneous Coulomb potential. In 1 + 1 dimensions we may represent the dirac 
matrices using pauli matrices, 7 = 03 and a = 7°7 1 = <j\. Despite the fact 
that Eq. (^) has no explicit lorentz covariance (space and time coordinates 
are treated differently in Eqs. (^,^)) the bound state energies for different 



cm. momenta k are correctly related: E = \/k 2 + M 2 , with M indepen- 
dent of k. In the limit of non-relativistic internal motion^ (e/m^ 1) the 
wave function xi x ) lorentz contracts in the standard way as a function of k. 
Related examples may be found in Ref. [[14] . 

The ground state wave function of QCD 4 is invariant under boosts. This 
is obviously not the case for our asymptotic states which according to Eq. 

3 This is in fact the only case where the solutions are normalizable and thus meaningful, 
due to the Klein paradox. The bound state momentum k can be arbitrarily large, however. 
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(p5[) contain bosons of definite 3-momenta (;$ Aqot). In the present formula- 
tion, the same boundary states must be used in all frames, since they model 
the same (invariant) ground state. Hence the perturbative expansion of a 
given QCD process will depend on the lorentz frame. If the method works, 
measurable quantities such as hadronic cross sections will be lorentz invari- 
ant. This does not mean that the hadron wave functions themselves will 
be invariant - the starting point of this paper was in fact that they appear 
phenomenologically to be strongly frame dependent. 

Acknowledgements. I would like to thank V. M. Braun and S. J. Brodsky 
for useful discussions. 
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